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Abstract. We show that all features of the ultrahigh energy cosmic ray spectrum from 
10^^ eV to 10^^ eV can be described with a simple power-like injection spectrum of 
protons under the assumption that the neutrino-nucleon cross-section is significantly 
enhanced at center of mass energies above ~ 100 TeV. In our scenario, the cosmogenic 
neutrinos produced during the propagation of protons through the cosmic microwave 
background initiate air showers in the atmosphere, just as the protons. The total air 
shower spectrum induced by protons and neutrinos shows excellent agreement with the 
observations. A particular possibility for a large neutrino-nucleon cross-section exists 
within the Standard Model through electroweak instanton-induced processes. 

1 Introduction 

The spectrum of cosmic rays extends in energy up to almost 10^^ eV. About 
twenty mysterious events were observed above 10^° eV by five different air 
shower observatories (AGASA [1], Fly's Eye [2], Haverah Park [3], HiRes [4], 
and Yakutsk [5]). Though some smaU-angle clustering in the arrival direction of 
the ultrahigh energy cosmic rays (UHECRs) is observed, the overall event distri- 
bution is isotropic. This indicates that they originate from several, isotropically 
distributed sources. 

Nucleons produced at large distances with energies above the Greisen-Zatse- 
pin-Kuzmin (GZK) cutoff [6] E'gzk- 4 • 10^^ eV interact with the cosmic mi- 
crowave background (CMB) and produce pions which decay into neutrinos. This 
way the nucleons lose their energy during propagation. The typical interaction 
length of nucleons above Eqzk is around 50 Mpc. Thus all events above 10^*^ eV 
should originate from small distances. However, no source within a distance of 
50 Mpc is known in the arrival directions of the post-GZK events. The angu- 
lar distribution of UHECRs above Eqzk does not show a correlation with our 
galactic plane which also indicates that they originate from large distances. No 
conventional explanation exists to the problem how can they reach us with en- 
ergies above 10^*^ eV without an apparent energy loss. 

At the relevant energies, among the known particles only neutrinos can prop- 
agate without significant energy loss from cosmological distances to us. It is this 

* Talk given at Beyond the Desert '03, Castle Ringberg, 9-14 June, 2003. 
** On leave from Institute for Theor. Physics Eotvos University, Budapest, Hungary. 



2 Z. Fodor et al. 



fact which led, on the one hand, to scenarios invoking hypothetical - beyond the 
Standard Model strong interactions of ultrahigh energy cosmic neutrinos [7] 
and, on the other hand, to the Z-burst scenario [8]. 

Interestingly, the flux of neutrinos coming from the pions produced during the 
propagation of nucleons - the cosmogenic neutrinos [7] - shows a nice agreement 
with the observed UHECR flux above Eqzk [9,10]. Assuming a large enough 
neutrino-nucleon cross-section at these high energies, these neutrinos could ini- 
tiate extensive air showers high up in the atmosphere, like hadrons, and explain 
the existence of the post-GZK events. This large cross-section is usually ensured 
by new types of TeV-scale interactions beyond the Standard Model, such as aris- 
ing through gluonic bound state leptons [11], TcV-scale grand unification with 
leptoquarks [12], or Kaluza-Klein modes from compactified extra dimensions [13] 
(see, however, Ref. [14]); for earlier and further proposals, see Refs. [15] and [16], 
respectively. 

In this review we discuss strongly interacting neutrino scenarios to solve the 
GZK problem, and in particular give an example which - in contrast to previous 
proposals - is based entirely on the Standard Model of particle physics. It exploits 
non-perturbative electroweak instanton-induced processes for the interaction of 
cosmogenic neutrinos with nucleons in the atmosphere, which may have a size- 
able cross-section above a threshold energy i?th = 0{{4:'Kmw / aw)^) / C^m-p) = 
0(10^^) eV, where mw denotes the W-boson mass and aw the electroweak fine 
structure constant [17-19]. We present a detailed statistical analysis of the agree- 
ment between observations and predictions from strongly interacting neutrino 
scenarios. 

Our scenario is based on a standard power-like primary spectrum of protons 

injected from sources at cosmological distances. After propagation, these protons 
will have energies below Eqzk, so they can well describe the low energy part of 
the UHECR spectrum. The cosmogenic neutrinos interact with the atmosphere 
and thus give a second component to the UHECR fiux, which describes the 
high energy part of the spectrum. The relative normalization of the proton and 
neutrino fluxes is fixed in this scenario, so the low and high energy parts of the 
spectrum are explained simultaneously without any extra normalization. Details 
of this analysis can be found in Ref. [20]. 

The structure of this review is as follows. In the next section we give the fluxes 
of protons and cosmogenic neutrinos both at their production and at detection. 
In Sect. 3 the possibility of using electroweak instantons as a source for large 
cross-section is discussed and the induced air shower rate is calculated. In Sect. 
4 we compare the predictions with observation and determine the goodness of 
flt, while conclusions are given in Sect. 5. 

2 Proton and cosmogenic neutrino fluxes 

We start with a power-like injection spectrum per co-moving volume of protons 
with energy Ei, spectral index a, and redshift (z) evolution index n, 

jp = jo Er" (1 + zf e{E^^^ - Ei) e{z - z^^) - z) . (1) 
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Here, jo is a normalization factor, Smax is the maximal energy, which can be 

reached through astrophysical accelerating processes in a bottom-up scenario, 
and Zmin/max takes into account the absence of nearby/very early sources. The 
overall normalization jo will be fixed by the observed flux, and our predictions 
are quite insensitive to the specific choice for i?max, z^im and 2;max, within their 
anticipated values. The main sensitivity arises from the spectral parameters a 
and n, for which we determine the 1- and 2-sigma confidence regions in Sect. 4. 

The propagation of particles can be described [9,21] by P^a{r, Ei; E) func- 
tions, which give the expected number of particles of type b above the threshold 
energy E if one particle of type a started at a distance r with energy Ei. With 
the help of these propagation functions, the differential Atix of protons (6 = p) 
and cosmogenic neutrinos (b = Vi^Vi) at earth can be given as 

F^iE) = l f dE, f dr (-)^^^l^ E,) . (2) 

In our analysis we go, according to dz = (1 + z)H{z)dr/c, out to distances 
-Rmax corresponding to ^max = 2 (cf. Ref. [22]), while we choose Zmin = 0.012 in 
order to take into account the fact that within 50 Mpc there are no astrophysical 
sources of UHECRs. We use the expression H^{z) — [Qm (1 + + i^A\ for 
the relation of the Hubble expansion rate at redshift z to the present one. Un- 
certainties of the latter, Ho ^ h 100 km/s/Mpc, with h = (0.71±0.07)xJ:^| [23], 
are included. Qm and Ha, with Qm + = 1, are the present matter and vac- 
uum energy densities in terms of the critical density. As default values we choose 
f^M = 0.3 and Ha = 0.7, as favored today. Our results turn out to be rather 
insensitive to the precise values of the cosmological parameters. 

We calculated Pb\a{r, Ei]E) in two steps, i) First, the SOPHIA Monte-Carlo 
program [24] was used for the simulation of photohadronic processes of protons 
with the CMB photons. For e+e" pair production we used the continuous energy 
loss approximation, since the inelasticity is very small (« 10~^). We calculated 
the P(,|a functions for "infinitesimal" steps (1 - 10 kpc) as a function of the 
redshift z. ii) We multiplied the corresponding infinitesimal probabilities starting 
at a distance r{z) down to earth with z = 0. 

Since the propagation functions are of universal usage, we decided to make 
the latest versions of —dPh\a/dE available for the public via the World-Wide- 
Web URL www.desy.de/~uhecr . 

3 Neutrino induced air shower rate 

The main assumption of our scenario is that the neutrino-nucleon cross-section 
al°^ suddenly becomes much larger than « 1 mb above center of mass energies 
^/s w 100 TeV. In this case, the corresponding neutrino interaction length = 
nip/al"^, with al°^ = al% + cr^^, falls below Xq = 1031 g/cm^ - the vertical 
depth of the atmosphere at sea level - above the threshold energy w 10^^ eV. 
Here and cr^jy denote the charged current and the new contribution to 
the cross-section. Above the neutrino threshold energy, the atmosphere becomes 
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opaque to cosmogenic neutrinos and most of them will end up as air showers. 
Quantitatively, this fact can be described by 



1 



■ A„(i3) 



(3) 



which gives the spectrum of neutrino- initiated air showers, for an incident cosmo- 
genic neutrino flux Fi, — '}2n[F,y- +-PpJ from Eq. (2), in terms of the atmospheric 
depth X{9), with 6 being the zenith angle. 
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Fig. 1. Left: Prediction of the electroweak instanton-induced neutrino- nucleon cross- 
section o"^^ (solid) in comparison with the charged current cross-section crj;'^ (dotted) 
from Ref. [25], as a function of the neutrino energy Ev in the nucleon's rest frame. Right: 
Neutrino interaction length due to combined effects of charged current interactions and 
instanton-induced processes. 



Such suddenly increasing cross-sections have been proposed in various mod- 
els involving physics beyond the Standard Model [11-16]. In Fig. 1 we show 
another example which is based entirely on the Standard Model exploiting non- 
perturbative electroweak instanton-induced processes [17-19]. Our quantitative 
analysis [20] was based on this cross-section, however it is expected to be insen- 
sitive to the exact form of it as long as it rises abruptly far above 1 mb. Note 
that such a behaviour is consistent with present upper bounds on electroweak 
instanton-induced cross-sections [26] . It should be noted that such a cross-section 
will lead, via dispersion relations, to lower energy deviations of Standard Model 
predictions for elastic scattering from their perturbative values. However, it is 
easily checked that, for the one shown in Fig. 1 (left), these corrections will be 
unobservably small in the energy regime available at present accelerators [27]. 

4 Comparison with UHECR data 

The predicted air shower rate induced by protons and neutrinos is given by 

^'prod(£'; a, n, £;,nax, ^min, -Zmax, jo) = Fp{E] . . .) + F^{E; . . .) . (4) 

We performed a statistical analysis to compare (4) with the observations and 
presented a measure for the goodness of the scenario [20] . We gave the best fit to 
the observations and the 1- and 2-sigma confidence regions in the {a,n) plane. 
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UHECR collaborations usually publish their results for the detected fluxes in 

a binned form. The first step of our analysis is to convert these fluxes into event 
numbers in each bin. We use the most recent results of the HiRes and AGASA 
collaborations and do our analysis separately with both data sets. We use the 
energy range 10^^'^ eV fO^^ eV which is divided into 38 equal logarithmic bins. 
In the low energy region, there are no published results available from AGASA 
and only low statistics results from HiRes-2. Therefore, we included the results 
of the predecessor collaborations - Akeno [28] and Fly's Eye, respectively - into 
the analysis. 

The goodness of the scenario is determined by a statistical analysis. We 
determined the compatibility of different (a.n) pairs with the experimental data. 
For some fixed {a,n) pair, the expected number of events in individual bins are 
(A = {Ai,...,Ar} with r being the total number of bins (in our case 38). The 
probability of getting an experimental outcome k = {ki,...kr} (where ki are 
non-negative integer numbers) is given by the probability distribution function 
P(k), which is just the product of Poisson distributions for the individual bins. 
It is easy to include also the « 30% overall energy uncertainty into the i-'(k) 
probability distribution. We denote the experimental result by s = {si,...Sr}j 
where the Sj-s are non-negative, integer numbers. The {a,n) pair is compatible 
with the experimental results if 



E 



P(k) < c . (5) 

k|P(k)>P(s) 

For a l-(or 2-)sigma compatibility one takes c=0.68 (or c=0.95), respectively. 
The best fit is found by minimizing the sum on the left hand side. 

Since we have 38 variables, it is practically impossible to calculate the sum 
in equation (5) exactly. If we rewrite it as 



E 



EkP(k)(?[P(k)-p(s) 



P(k) = ypV) ' (6) 

k|P(k)>P(s) ^ ' 

then we can calculate this sum approximately using an importance sampUng 
based Monte-Carlo technique. We have to generate the components of k with 
Poisson distribution and take only those in the sum, for which -P(k) > P(s). 

Figure 2 (left) shows our best fits for the AGASA and for the HiRes UHECR 
data. The best fit values are a = 2.68(2.68) and n = 2.65(2.9), for AGASA(HiRes). 
We can see very nice agreement with the data within an energy range of nearly 
four orders of magnitude. The fits are insensitive to the value of .Emax as far as 
we choose a value above ~ 3 • 10^^ eV. The shape of the curve between 10^^ eV 
and 10''^^ eV is mainly determined by the redshift evolution index n. At these 
energies the universe is already transparent for protons created at z « while 
protons from sources with larger redshift accumulate in this region. The peak 
around 4 • 10^" eV shows the accumulation of particles due to the GZK effect. 
Neutrinos start to dominate over protons at around 10^° eV. 

It is important to note that, if we omit the neutrino component, then the 
model is ruled out on the 3-sigma level for both experiments. This is due to the 
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Fig. 2. Left: Ultrahigh energy cosmic ray data (Akcno + AG AS A on the upper panel 
and Fly's Eye + HiRes on the lower panel) and their best fits within the electroweak 
instanton scenario (solid) for iSmax = 3 • 10^^ eV, Zmin = 0.012, Zmax = 2 consisting of a 
proton component (dotted) plus a cosmogenic neutrino-initiated component (dashed). 
Right: Confidence regions in the a-n plane for fits to the Akeno + AGASA data (2- 
sigma (long dashed)) and to the Fly's Eye + HiRes data (1-sigma (solid); 2-sigma 
(short-dashed)), respectively. 



fact that we excluded nearby sources by setting z^in 7^ (see also Ref. [29]). A 
choice of Zmin = makes the HiRes data compatible with a proton-only scenario 
on the 2-sigma level (see also Rcfs. [4,30]). 

Figure 2 (right) displays the confidence regions in the {a,n) plane for AGASA 
and HiRes. The scenario is consistent on the 2-sigma level with both experiments. 
For HiRes, the compatibility is even true on the 1-sigma level. It is important 
to note that both experiments favor the same values for a and n, demonstrating 
their mutual compatibility on the 2-sigma level (see also Ref. [31]). 

Finally, let us discuss the consistency of our scenario with the currently avail- 
able limits on deeply penetrating showers from Fly's Eye [32] and AGASA [33]. 
Taking into account - in distinction to Ref. [34] - the atmospheric attenuation 
of the cosmogenic neutrino flux predicted in our scenario and the uncertainties 
in the estimate of the range of depth within which the shower must originate to 
trigger the array, we find that AGASA should have seen 1-10 quasi- horizontal 
air showers (^^60°) from the electroweak instanton-induced processes during a 
running time of 1710.5 days. This is consistent with AGASA's present analysis 
of their respective data [33] . The Fly's Eye upper limit on the product of the to- 
tal neutrino flux times neutrino-nucleon cross-section, (F,y a^J'j^) Fly's Eye [32], in 
the energy range 10^^ ~ eV, can be translated, for a given predicted neutrino 
flux FP'-<=d, into an upper hmit on al% < (F^ cr^,°7^)Fiy's Eye/i^P"='^, as long as it 
is smaller than 10 [18,35]. We find that, for our predicted cosmogenic neu- 
trino flux, the right-hand-side of this inequality is larger than 10 fib in the whole 
energy range, such that the Fly's Eye non-observation of quasi-horizontal air 
showers does not give any constraint. We therefore conclude that our prediction 
of the neutrino-nucleon cross-section, as shown in Fig. 1 (left), does not contra- 
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diet any constraints from cosmic ray experiments so far, as long as the ultrahigh 
energy cosmic neutrino flux is at the cosmogenic level we have predicted. 

5 Summary and conclusions 

We have shown that a simple scenario with a single power-like injection spectrum 
of protons can describe all the features of the UHECR spectrum in the energy 
range 10^^ ~ eV. In our scenario, the injected protons produce neutrinos dur- 
ing their propagation and these neutrinos arc assumed to have large enough 
cross-section to produce air showers high up in the atmosphere. As an example 
we discussed the possibility that Standard Model electroweak instanton-induced 
processes may give a cross section which is suitable for this scenario. The model 
has few parameters from which only two - the power index a and the redshift 
evolution index n - has a strong eff'ect on the final shape of the spectrum. We 
found that for certain values of a and n this scenario is compatible with the 
available observational data from the AGASA and HiRes experiments (com- 
bined with their predecessor experiments. Fly's Eye and Akeno, respectively) on 
the 2-sigma level (also 1-sigma for HiRcs). The ultrahigh energy neutrino com- 
ponent can be experimentally tested by studying the zenith angle dependence 
of the events in the range 10^^ ~ eV and possible correlations with distant 
astrophysical sources [36] at cosmic ray facilities such as the Pierre Auger Ob- 
servatory [37], and by looking for enhanced rates for throughgoing muons at 
neutrino telescopes such as AMANDA [38] . As laboratory tests, one may search 
for a model-independent enhancement in (quasi-)elastic lepton-nuclcon scatter- 
ing [27] or for signatures of QCD instanton-induced processes in deep-inelastic 
scattering [39], e.g. at HERA. 
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